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We propose a simple mechanical model of locomotion induced by internal force and friction. We 
first construct a system of two elements as an analog of the bipedal motion. The internal force 
does not induce a directional motion by itself because of the action-reaction law, but a directional 
motion becomes possible by the control of the frictional force. The efficiency of these model systems 
is studied using an analogy to the heat engine. As a modified version of the two-elements model, we 
construct a model which exhibits a bipedal motion similar to kinesin's motion of molecular motor. 
Next, we propose a linear chain model and a ladder model as an extension of the original two-element 
model,. We find a transition from a straight to a snake-like motion in a ladder model by changing 
the strength of the internal force. 
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^ \ PACS numbers: 05.65.-|-b, 45.05.-l-x, 87.19.lu 

Q , I. INTRODUCTION 

D 

A uni-directional motion is one of recent topics in nonequilibrium physics. The directional motion has been in- 
' tensively studied, motivated by molecular motors [H-Q- Molecular motors such as kinesin and myosin exhibit uni- 
, directional motion along microtubules and actin filaments Q. Molecular motors are interpreted as energy converters 
C/3 from chemical energy to mechanical energy. A ratchet mechanism is considered to be a possible mechanism to generate 
^ , a directional motion using thermal energy or chemical energy [gI, (tJ . 

' Directional motions have been studied also in physical and chemical systems. Linke et al. demonstrated experi- 
mentally that liquid drops perform a self-propelled motion when they are placed in contact with hot surfaces [8|. A 
sawtooth type asymmetric surface played an important role for the directional motion in their experiments. Chaud- 
hury and Whitesides reported an uphill motion of a liquid drop induced by a surface chemical gradient @. Sumino 
Q , et al. found a reciprocal motion of an oil drop on a glass plate induced by oscillatory chemical change of the surface 
' tension [lo| . A directional motion was also observed in the active gels where the mechanical motion is coupled with 
oscillatory chemical reactions (ill. IT^. In the experiment by Maeda et al., the time-periodic volume change induced 
by the chemical reaction is transformed to a directional motion on an asymmetric surface. We studied a directional 
motion in a simple model of a self-oscillating elastic filament, motivated by the experiments of the active gels [TS]. 
■ Toyobe et al. performed an experiment of a microscopic particle climbing upwards on a spiral-stair-like potential 
CD ' using a feedback control 

04 . The self-driving mechanism, especially, the bipedal motion has been one of the main topics in robotics 15]. The 
well-controlled walking robots have been developed, where the positions and motions of legs, feet, and many joints are 
, precisely controlled by micro computers. However, it is known that there is another simple walking mechanism called 
^D ' passive dynamic walking, where no intricate control is needed jT6| . Bipedal walking is realized along a downward 
slope by a simple mechanical mechanism. As for the locomotion of animals, there are a large amount of literatures 
in biology. For example, the reptation motion of snakes was mathematically discussed in [17|. The energetic cost 
of the snake motion was estimated [isj . Recently, snake- type robots were constructed by combining biology and 
robotics 

^ In this paper, we propose simple mathematical models based on elementary mechanics. Directional motion is 

, induced by periodic change of internal force and friction in our model. We study the simple mechanical model 
systems from a view point of the nonlinear-nonequilibrium physics, hoping that the simple models might be useful 
to understand directional motions ranging from molecular motors to animal motion and robotics. In §2, we study 
first a simple system composed of two elements, demonstrate a directional motion and discuss the efficiency as an 
engine. Then, we modify the two-element model to apply to kinesin's motion, and propose Langevin type equations, 
because thermal noises are expected to play an important role in molecular motors. In §3, we extend the original 
model studied in §2 to a linear chain model. In §4, we generalized the linear chain model to a ladder model where two 
linear chains are coupled, and we find a dynamical instability from which a snake motion appears. The snake motion 
is usually considered to be generated by internal muscle forces but our model system suggests that a snake motion 
might be generated from a dynamical instability, when the body size is sufficiently contracted during a uni-directional 
motion. We consider a sliding friction model and a viscous friction model respectively for a bipedal model and a linear 
chain model. The sliding friction model might be more realistic for the directional motion on the ground. However, 
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some analytical predictions are possible for the viscous friction models, since the simple system is a linear system, 
and the viscous friction models might be useful for theoretical understanding of the uni-directional motion. 

II. BIPEDAL MOTION IN A SYSTEM OF TWO ELEMENTS 

We first study a system composed of two elements as an analogy of the bipedal walking motion. In the bipedal 
motion, two legs move and stop alternatively. The motion is driven by some internal muscles and stops by the friction 
with the ground. The frictional force is proportional to the normal force on the ground. The normal force changes 
periodically by the positional shift of the center of mass of the body. That is, one leg moves forward when the weight 
of the body is loaded on the other foot, and the roles of two legs are reversed in the next half period. The actual 
bipedal motion is rather complicated, including intricate controls of several joints and muscles. Here, we consider a 
very simple mathematical model based on the elementary mechanics: 

K{x2 — xi — a + Fsinut) — W — Fri, 

-K{x2 - xi - a + Fsinuit) -W - Fr2, (1) 

where xi,2 is the position of the two elements (analogues of legs), m is the mass, K denotes the spring constant of a 
spring between the two elements which is an analogue of muscles, a is an average interval between the two elements, 
and the natural length of the spring is assumed to change in time as a — Fsinuit by the internal force. W is the 
external load. In a climbing motion on a slope of angle 6, the external load is represented as = m^sin^. Fri, Fr2 
are the frictional forces. A typical time scale in this model system is y^m/K, which is proportional to the period of 
the oscillation of the spring. If the time is rescaled with the typical time ^/m/K, and other parameters are rescaled 
such as W/K W, Fri/K Fr^,Fr2/K Fr^ and {m/Ky/'^uj to, Eq. (1) is rewritten as 

{x2 — a;i — a + F sin uit) — W — Fri , 

-{x2 - xi - a + Fsincot) -W - Fr2, (2) 

If the friction and the external load are absent, the center of mass {xi + X2)/2 does not move if the initial velocities 
are zero, because of the action-reaction law. We further assume the sliding friction Fri and Fr2 of the form: 

Fri = - /3sina;t)sgn(t'i), 

Fr2 = Ai<i(l + /3sinwi)sgn(w2), (3) 

where /x^ is the coefficient of the sliding friction, and /3 is the amplitude of the periodic variation of the normal force, 
and Vi = dxi/dt,V2 = dx2/dt. If the motion stops, the static friction works and the maximum static friction force is 
expressed as /is(l T (Ssinuit). Note that the time- variation of the normal force for the two elements is out of phase, 
which comes from the shift of the center of mass of the body. The average interval a between the two elements is 
irrelevant in this model. The system of a = might be better for the bipedal walking, and the model with positive a 
might correspond to a system composed of fore and hind legs. 

At first we show a numerical result for (3 = 0, that is, the normal force is constant in time. Figure 1(a) shows the 
time evolution of Xi{t) and X2{t) for a = 1.5, F = 0.5, W = 0,u) = 27r/60,^ = 0,1x4 = 0.25 and Hs = 0.6. Because 
of the periodic change of the internal force, the interval of X2 — xi exhibits oscillation. A stick-slip motion of stretch 
and contraction appears owing to the alternation of the sliding and static frictions. However, the center of the mass 
(a;i +X2)/2 does not move. For (3 = 0.8, a uni-directional motion appears as shown in Fig. 1(b). Sliding motion 
occurs alternatively for Xi{t) and X2{t). When the spring is contracted, the normal force on the foreleg is sufficiently 
large and the foreleg cannot move owing to the strong static friction, however, the normal force on the hind leg is 
relatively weak, and the hind leg can slide forward. On the other hand, when the spring is stretched after a half period, 
the normal force on the foreleg is weak and the foreleg slides forward, however, the normal force on the hind leg is 
sufficiently large, and the hind leg cannot move. That is, the hind leg moves forward fixing the position of the foreleg 
in the first half period, and the foreleg moves forward fixing the position of the hind leg in the second half period. 
The internal forces work in the opposite directions for the fore and hind legs because of the action-reaction law of the 
internal force, but the directional motion is induced owing to the timing between the internal force and the temporal 
change of the friction force. This is a physical mechanism of the bipedal motion in our model. Even if the external 
load W is not zero, a uni-directional motion in the x-direction is induced, which corresponds to the bipedal walking 
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FIG. 1: (a) Time evolution of xi{t) and X2{t) for Eqs. (2) and (3) at a = 1.5, F = 0.5, W ^0,uj = 27r/60, ^ = 0, /id = 0.25, and 
Us = 0.6. (b) Time evolution of xi(t) and X2{t) at a = 1.5, F = 0.5, W ^ 0,lj = 27r/60, /3 = 0.8, fid = 0.25, and = 0.6. (c) 
Ji = {dxi/dt + dx2/dt) vs. W at a = 1.5, F = 0.5, uj = 27r/60, /? = 0.8, /id = 0.25 and fj.s = 0.6. 



along an upward slope. Figure 1(c) shows a long time average of the sum of the two velocities: Ji = {dxi / dt + dx2 / dt) 
as a function of W. The other parameter values are set to be a = 1.5, F — 0.5, u — 27r/60,/3 — 0.8, /id = 0.25 and 
/is = 0.6. The average velocity decreases with the external load W. There is a discontinuous transition to the state 
of Ji = at = 0.45 in this model. The uni-directional motion is impossible for W > 0.45. The discontinuous 
transition occurs because the frictional force is a nonlinear function of the velocity. 

The dynamics is not easily treated mathematically in the nonlinear system. We can propose a linear model with 
viscous friction: 
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{x2 — xi — a ^- F sin ujt) — W~Fri, 
— (x2 — Xi — a + Fsmut) — W — Fr- 



dt^ 

Fri — /i(l — /3 sina;i)wi, 

Fr2 = /i(l + /3sina;i)u2, (4) 

where /i is the coefficient of the viscous friction. Figure 2(a) is the time evolution of xiit) and X2{t) for a = 1.5, F = 
1,W = Q,ljJ — 27r/60,/3 = 0.5, and /i — 20. A forward motion is observed on the average as in Fig. 1(b) even in this 
simpler model, however, the stick-slip motion does not appear in this model. Although the internal force works in the 
opposite directions for the two elements, a directional motion appears because of the time-periodic frictional force. 
This might be interpreted as a kind of ratchet mechanism, in that the reverse motion is effectively suppressed. The 
direction of the motion is determined by the timing of the variation of the internal force and the frictional force. The 
energy is injected through the forcing term of F sinujt, and is dissipated by the frictional force in this system. That 
is, our model system is a typical nonequilibrium system. In a molecular motor like the kinesin, the ATP is a source 
of energy and the energy is dissipated most through the viscous friction with water molecules. 

Since the model equation (4) is a simple linear system, an approximate solution can be obtained by assuming 

xi{t) — v^t + Aainut -\- B cosujt, 

X2{t) — a + VQt — A sin Ldt ~ B cos ujt. (5) 

Substitution of Eq. (5) into Eq. (4) yields 

-Lj'^A = -2A + /.iLoB + F + fil3vo, 
-Lo'^B = -2B - fiojA, 

vo = -W/fi- PujB/2, (6) 

if the terms including sin(2a;t) and cos(2a;i) are neglected. The higher harmonics appear through the terms of 
— /i/3(sina;i)z; in Fri and Fr2, and such higher harmonics can be neglected when /3 is sufficiently small. The solution 
of Eq. (6) is 

^_ {u^-2){pW-F) 



(a;2-2)2 + /i2w2(i_^2/2)' 
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FIG. 2: (a) Time evolution of xi{t) and X2{t) for Eq. (4) at a = 1.5, F ^ 1,W ^ 0,uj ^ 27r/60, /? = 0.5 and fj, = 20. The dashed 
curves denote Eq. (5), whose parameters are given by Eq. (7). (b) L12 (solid line) and L21 (dashed line) for the linear model 
(4) as a function of /3 for a = 1.5, uj = 27r/60 and \i = 20. The solid and dashed lines are completely overlapped. The dotted 
line denotes L12 = L21 by Eq. (9). (c) Efficiency 77 (solid line) at the maximum power for a = 1.5, uj — 27r/60 and /i = 20. The 
dashed line is the approximation by Eq. (9). (d) Efficiency rj (solid line) at the maximum power for a = 1.5,aj = 27r/20 and 
fj, = 20. The dashed line is the approximation by Eq. (9). 



(a;2-2)2 + /i2^2(i_^2/2)' 

"° = -^/^ 2{(C.2- 2)2+^2^2(1,^2/2)} 

A divergence occurs at uj = \/2 and /3 = in Eq. (7). The frequency w — \/2 is the resonance frequency of the 
Hnear vibration. In our numerical simulation, uj < \/2 is assumed, i.e., the forcing period is sufficiently longer than 
the period of the proper vibration, and /? < 1 is assumed because the normal force becomes negative if /3 > 1. 
Therefore, the divergence does not occur in our parameter setting. For F = 1,W ~ 0,uj = 27r/60, (3 = 0.5 and ^ = 20, 
A = 0.255, B = —0.269 and ?;o = 0.00703 are obtained. The dashed curves in Fig. 2(a) represent the approximate 
solutions (5). Fairly good agreement is seen. 

As suggested from Eq. (7), a uni-directional motion in the x-direction is induced also in this viscous friction model, 
even if W is slightly positive. When W is positive, the uphill occurs, which implies that this system can perform 
the "work" in the classical mechanics. If this system is interpreted as an energy converter from mechanical energy 
to potential energy Vl^(a;i + X2), the theory of the heat engine can be applied to this system [g^, [2l|- The efficiency 
of the heat engine at the maximum power is recently discussed in several thermodynamic systems. The efficiency at 
the maximum power is practically useful because the Carnot's maximum efficiency is attained only in the quasi-static 
process Van den Broeck applied Onsager's reciprocal relation to the efficiency of a heat engine at the maximum 
power. The concept of the efficiency at the maximum power has not yet been studied in detail in deterministic 
systems. In Onsager's theory, two thermodynamic flows Ji and J2 are expressed with the linear functions of two 
thermodynamic forces W and F. In our system, two kinds of flows Ji = vi + V2 and J2 = {(vi — V2) sinwt) can be 
deflned, where (• • •) represents a long-time average. Ji represents the flow of mass. The input energy by the internal 
force per unit time is written as -B = {viF sinujt + V2{— F sin ujt)) ~ FJ2. J2 represents a thermodynamic flow with 
respect to the input energy. In the above approximation of Eq. (5), Ji — 2vq and J2 = —Buj are obtained, and then, 
the Onsager type relations are derived in this deterministic system as 

Ji = Lii{-W) + L12F, 

J2 = L2i{-W) + L22F, (8) 

where 
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FIG. 3: (a) Time evolution of xi{t) (solid curve) and X2{t) (dashed curve) for F=la,tg = 5,l3 = 0.9, — 2-k/6Q,K = 2 
and = 10 in Eq. (10). (b) Time evolution of xi{t) (solid curve) and X2(t) (dashed curve) for F = 2. (c) Average velocity v 
as a function of F at g = 5, ^ = 0.9, lj = 27r/60, K = 2 and fi = 10. (d) Time evolution of X2{t) for F = 0.2 and T = 3 at 
g ^5,13 ^ 0.9, cj = 27r/60, isT = 2 and = 10. 



In this approximation, Onsager's reciprocal relation L12 = L21 is satisfied. The power, that is the work per unit 
time, is expressed as P = WJi = —LuW^ + L12WF, and the input energy by the internal force is written as 
E — FJ2. The dissipated energy by the friction or the waste heat per unit time is expressed as Q = E — P = 
Lii{-Wf + (L12 + L2i){-W)F + If L21 = L12, Q = Lii{-Wf + 2Li2{-W)F + L22P^ > 0. The efficiency 

rj is expressed as ry = P/ E, which is the ratio of the work per unit time over the input energy per unit time. The 
maximum power is obtained when dP/dW — 0, that is, W = {L12F) / {2Lii). The efficiency at the maximum power 
is expressed as 77 = Lf2/(4L22iii — ^12^21)- 

Figure 2(b) shows numerically estimated values of L12 (solid line) and L21 (dashed line) for the linear model (4) as 
a function of /3 at a = 1.5, a; = 27r/60 and /i — 20. The dotted line denotes L12 — L21 in Eq. (9). The approximation 
is rather good for /3 < 0.7 and Onsager's reciprocal relation is satisfied in the whole range. Figure 2(c) shows the 
efficiency 77 at the maximum power. The dashed line is 77 = L^2/(4-^22iii — ^12^21) obtained by using Eq. (9). Good 
agreement is seen for the efficiency for (3 < 0.7. For /3 > 0.7, higher harmonics become important, and the deviation 
becomes large. The efficiency depends on ^ and w, i.e., rj increases with fi and decreases with uj. Figure 2(d) shows 
the efficiency 77 at the maximum power as a function of /? for a = 1.5, uj ~ 27r/20 and [i — 20. The period l-n juj is one 
third of the case of Fig. 2(c) and the efficiency is fairly increased. 

We can modify the previous model to apply to bipedal motion in molecular motors. It is known that the kinesin 
exhibits a bipedal motion along the microtuble. It is characteristics of the system that the microtuble has a spatially- 
periodic structure of period Iq. The one leg of the kinesin becomes docked to the microtuble owing to the bound ATP. 
The other leg moves forward by 2^o over the docked leg. At the next half period, the other leg becomes docked and 
the former leg moves forward by 11^ over the other docked leg. Stepwise motions occur alternatively for the two legs. 
Thermal fluctuations play an important role in molecular motors because of the small size. The effect of thermal 
noises can be taken in by Langcvin ype equations. We propose a model for a bipedal motion of the molecular motor 
along a spatially periodic potential as 

K(x2 — xi + Fsinwt) — ^ fismuji) sin27ra;i + ^i(t), 

(ix2 

—K(x2 — Xi + Fsincjf) — /i 1" 5(1 + sinwi) sin 27ra;2 + ^2(0, (10) 

di 

where the spatial period is assumed to be 1. ^i(t) is assumed to be Gaussian white noises which satisfy {^i{t)^j {t')) ~ 
2TSi,j6{t — t') for i,j — 1 or 2. The parameter T is proportional to the temperature. If T = 0, Eq. (10) is reduced to 
be a deterministic equation. The parameter a is set to be 0, because there is no difference of the fore and hind legs 
in the kinesin system. The amplitude of the spatially-periodic potential cos 27ra;/(27r) is assumed to change in time 
and the timing is out of phase for the two elements. The fact that the attractive force between the kinesin and the 
microtuble changes in time by the binding of ATP is simply represented by the time-periodic function 1 =F /5 sin ujt in 
our model. The two elements tends to be docked near potential minima at n + l/2 where n is an integer. On the other 
hands, the internal force Fsinut tends to enforce exchange of position for the two elements. We have first performed 
numerical simulation of Eq. (10) for T = at 5 = 5, /? = 0.9, uj = 27r/60, K = 2, and /i = 10. Figure 3(a) displays time 
evolutions of xi{t) (solid curve) and X2{t) (dashed curve) for F = 1. Stepwise motion appears for xi{t) and X2{t). 
When xi stays near a potential minimum, X2 moves forward by 21q = 2. On the other hand, when X2 stays near a 
potential minimum, xi moves forward by 21q — 2. The stepwise motion occurs alternately for xi and X2- The average 
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FIG. 5: (a) Time evolution of the profile Vi iox N = 2h,W = F = l,k = lTY/N,e = 7r/2,a = l,uj = 27r/50,/i = 15 
and /? = 0.2 in the viscous friction model (12). (b) Onsager's coefficients L12 (+) and L21 (rhombus) as a function of /3 for 
A'^ = 25, fe = 4:n/N,6 — 7r/2,cj — 2-k/50 and /i — 15. The dashed lines are the approximation by Eq. (14). (c) Efficiency 
r; as a function of /? by numerical simulation and the approximation using Eq. (14). (d) Efficiency ?7 as a function of 6 for 
A'' = 25, fc — A-k/N,uj — 27r/50,pi = 15 and /3 — 0.6 by numerical simulation and the approximation using Eq. (14). 



velocity of the stepped motion is 2lo/{2n/uj) = 2/60 = 0.0333. Thus, we have succeeded in reproducing qualitatively 
a bipedal motion of kinesin. When F is smaller than 0.77, such a stepped motion disappears and the two elements are 
oscillating around a potential minimum. When F is larger than 1.94, a stepped motion with stepping length AIq = 4 
appears as shown in Fig. 3(b) for F — 2. The average velocity is 0.0666 for this stepped motion. Figure 3(c) shows 
the average velocity v a,s a, function of F. The average velocity changes stepwise. It is a phenomenon analogous to 
the frequency locking. 

The effect of thermal noises was studied at g = 5,/3 = 0.9, w = 27r/60,ii' = = 10,F 0.2 and T = 3. The 
stepwise forward motion does not occur and the average velocity is zero at F = 0.2 for T = as shown in Fig. 3(c). 
Figure 3(d) shows time evolution of X2{t) at T = 3. A noisy but stepwise forward motion appears in Fig. 3(d). One 
step is = 1 in most cases. It implies that the forward motion is induced by the thermal noises. Our model is similar 
to thermal ratchet models studied in [l]-[4], especially similar to a flashing ratchet model ^], however, our model is 
unique in that the potential is not a sawtooth type but mirror-symmetric. The attractive force to the spatially-periodic 
filament changes in time and it is out of phase for the two elements, which is essential for the directional motion. 

III. LOCOMOTIVE MOTION IN A LINEAR CHAIN MODEL 

In this section, we consider a one-dimensional chain model as a simple model of an earthworm-type motion. It is 
an extension of the two-element models expressed by Eq. (2) with Eq. (3), and Eq. (4). A sliding friction model is 
written as 

— {xi+i — Xi — a + Fsm(iut + a^)} — {xi — Xi^i — a + F sm( Lot + a,;_i)} — Fri — W, 
Fri = - /?sin(wt + ai + 6')}sgn(ui), (11) 
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where i = l,2,---,iV, A'^ is the total number of elements, ai = ki, and a parameter 9 represents a phase shift 
between the internal force and the normal force. In this model, the internal force and the normal force are set to 
propagate like a wave along the linear chain. The parameters w and k denote the frequency and the wavenumber of 
the wave. For an earthworm, such a wavy deformation of body segments is actually observed. The normal force on 
the ground is temporally controlled by the local change of the angle of tiny bristles on the skin. The configuration 
of the normal force along the long body can be widely changed by a slight change of the anchoring state to the 
ground, because the normal force tends to be concentrated near locally bristled regions. We have further assumed 
the free boundary conditions at z = 1 and i = N. Figure 4(a) displays time evolution of Xi{t) for z = 10 at 
iV = 25,a = 4:'K/N,e tt/2,F = 0.5, a; = 27r/50,/3 = 0.2, fid = 0.25 and ^ 0.6. A stick-slip motion is 

observed owing to the alternation of the sliding and the static friction. Figure 4(b) displays time evolution of profile 
Vi for the same parameter set. The stick-slip motion propagates from the head {i = N) toward the tail {i = 1) like an 
earthworm. The whole body moves forward successively by the propagation of the locally sliding region from head to 
tail. 

The viscous friction model is written as 

= {xi+i — Xi ~ a + Fsm{u!t + ai)} — {xi — Xi-i — a + Fsin(a;t + ai-i)} — Fri — W, 



Fri = IJ.{1- |3sin{uJt + a^ + e)}vi. (12) 

Since this viscous friction model is a linear system, it can be solved approximately. If Xi{t) = VqI + Asmiut + ai) + 
Bcos{ujt + ai) is assumed and the higher harmonics are neglected, the amplitudes A, B and the average velocity vo 
are given by 

{liJ-ic siu cori 0/2 -iv- + 2-2 cos A- jFi + fiu.-( l - f cos- 61/2)^2 



B = 



(w2 + 2 cos fc - 2)2 + ^2^2 _ ^2^2^2/2 

-^lUJ{l - /32 sin2 e/2.)Fi + {-plP^lo sin 6* cos 6*72 - u;2 + 2 - 2 cos fcjfa 



(u;2 + 2 cos fc - 2)2 + ^2^2 _ ^2^2^2/2 

Vo = I3uj AsmO /2 - liujB cose /2-W/n. (13) 

where Fi = (1 — cos k)F — Wf3 cosd and F2 = Fsink — W/SsinO. Figure 5(a) shows time evolution of the profile 
of Vi{t) N = 25,W = 0,F = 1,9 = 7r/2,fc = 47r/7V,a = l,w = 27r/50,Ai = 15 and /3 = 0.2. The deformation 
propagates smoothly from the head {i = N) toward the tail {i = 1) and the stick-slip motion does not appear in 
this model in contrast to the sliding friction model. The solid curves denote numerical results and the dashed curves 
express the theoretical result using Eq. (13). Fairly good agreement is seen except near the two boundaries i = 1 and 
N. 

We can discuss the efficiency also in this model. Similarly to the two-element model, two types of flows are 
defined as Ji = ^-"^^ "^2 — ^ Wi+i)sin(wt + 0;^)). The two types of flows are expanded as 

Ji = N{L-i_i{-W) + L12F} and J2 = N{L2ii-W) + L22F}. The coefHcients Ln ~ L22 are expressed by Eq. (13) as 



ill 
L12 
L21 



/i 2{ (a;2 + 2 cos A: -2)2 + ^2^2 _ ^2^2^2/2}' 

/3w[(— w2 -)- 2(1 — cosfc)){(l — cos A:) sin^ — sin fc cos 6*} + w/x{(l — cosfc) cos^ + sin A; sin ^}] 

2{ (a;2 + 2 cos fc - 2)2 + /x^w^ - f^^p^cv^ /2} 
/3a;[(— ^2 -|- 2(1 — cos A;)){sinfccos6' — (1 — cos k) sin 9} + cj/i{(l — cos A;) cos^ + sin A; sin ^}] 

2{(cj2 + 2 COS A: - 2)2 + /i2cj2 _ ^2^2/2} ■' 



_ LU^ lUslsm"^ k + {1 — cos k)^} — u)'^nP'^/2{siiik cos 9 — {I — cos k)siii9}'^ 

" 2{(w2 + 2cos(A;) - 2)2 + ^I'^iv^ - /i2^2^2/2} ^ ' 

Figure 5(b) shows L12 (+) and ^21 (rhombus) as a function of (3 iov N = 25, k = 4:Tr/N,9 = tt/2,ui = 27r/50 and 
/i = 15. The dashed and dotted curves are respectively L12 and L21 by the approximation (14). The approximation 
is good for small (3, but the deviation becomes large for large /3 owing to the higher harmonics. Onsager's reciprocal 
relation L12 = L21 is not exactly satisfied in this system. It may be related to a fact that the time-reversal symmetry 
is not satisfied in the model equation (12). We do not understand well, but Onsager's reciprocal relation might be 
satisfied in the bipedal model (4) because of an additional symmetry with respect to the exchange of xi and X2- The 
efficiency r] at the maximum power is evaluated as rj = L\2/ {^L22L\\ — L12L21) even in the case of ii2 ^ ^21 in this 
linear system. Figure 5(c) displays numerically obtained values of 77 as a function of /3 and the approximate values using 
Eq. (14). The efficiency 77 depends on other parameters. Figure 5(d) displays numerically obtained values of 77 as a 
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(o) (d) 




FIG. 6: (a) Schematic figure of a ladder model, (b) Five snapshot patterns of the ladder model for A'" = 50, F — 0.8, a = l,k — 
47r/25, 8 — 7v/2, uj — 27r/50, j3 = 0.9 and ^ = 15 in the viscous friction model (15) and (16). (c) Mean square (Sy^) as a function 
of F. (d) Average body length L as a function of F. 



function of the phase difference 9 and the approximate values using Eq. (14) for A'' = 25, k — Att/N, uj = 27r/50, /i = 15 
and /3 — 0.6. The efficiency takes a maximum near 9 = 7r/2, and almost near 9 = and tt. This figure implies 
that the timing between the internal force and the frictional force is very important for the locomotive motion in this 
model. 



IV. REPTATION MOTION IN A LADDER MODEL 



We can generalize the linear chain model in the previous section to a ladder model where two chain models are 
laterally coupled as shown in Fig. 6(a) as a very simple model of a snakelike motion. In the snakelike motion, the 
body wriggles in the y-direction, when it moves in the a;-direction on the average. The flex deformation appears along 
the body in the snakelike motion. A linear chain exhibits no resistance against the flex deformation. In a ladder 
model, the fiexural rigidity becomes nonzero owing to the finite thickness. This is a reason why we have introduced 
the ladder model. 

Although the snaking motion was observed in a sliding friction model, we show numerical results of a viscous friction 
model for the sake of simplicity. The viscous friction model is written as 

= ^(|xj-x.,|-aj,)ej, -^Vi, (15) 
j 

F„ = l3 sin{ujt + ai + 9)}v,. (16) 

where = {xi,yi) is the two-dimensional position of the ith element, ai = ki and a parameter 9 denotes a phase 
difference between the frictional force and the internal force. The elements of i = 1, 2, • • • , A^/2 are located in the 
lower chain, and the elements of i = A^/2 + 1, • • • , A^ are located in the upper chain. The summation in Eq. (15) 
is taken for the five neighboring sites of i, eji = (x^ — Xi)/|xj — Xi|,aji = a — F sin(ujt + ai) for the right site, 
ttji = a — Fsin(aji + ai-i) for the left site, aji = a for the opposite site, aji = yja? + {a — Fsin(a;i + ai)}^ for the 
diagonally right site, and a^i = ^Ja^ + {a — Fsin(a;t + a^-i)}^ for the diagonally left site. Figure 6(b) shows five 
snapshot patterns of the ladder model for A^ = 50, F = 0.8, a = l,fc = 47r/25,a; = 27r/50,/3 = 0.9 and /it = 15. 
The time interval between the two snapshots is 800. In this simulation, the y coordinate yN/2ii) and yN{t) at the 
right boundary are respectively fixed to 10 and 11. This boundary conditions induce a uni-directional motion in the 
x-direction on the average. If the free boundary conditions are imposed at the right boundary, the direction of motion 
changes in time. The ladder moves forward in the x-direction and it is deformed in the y-direction. As a result, 
a snakelike motion appears. We have calculated the mean square {Sy"^) = {yf) — {yi)"^ to quantify the deformation 
in the y-direction. The average with respect to i and the long-time average are taken for the calculation of {Sy"^). 
Figure 6(c) shows {5y'^) as a function of F . The mean square {Sy"^) changes continuously from to nonzero values 
a.t F = Fc ^ 0.52, which implies the appearance of the snaking pattern. The ladder is straight and exhibits an 
earthwormlike linear motion for F < Fc. The straight state becomes unstable and the snake pattern appears for 
F > Fc- Figure 6(d) shows the average body length L = + — Xi — xj^/2+i)/2 as a function of F. The 
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FIG. 7: (a) A snapshot pattern of the ladder model for iV = 150, F = 0.69, a = l,k = 47r/25, 6 ^■k/2,u = 27r/50, ^ = 0.9 and 
= 15 in the viscous friction model, (b) A snapshot pattern of the ladder model for k = 27r/25. The other parameters are the 
same as (a), (c) Wavenumber q of the snaking pattern as a function of k for F = 0.4 and = 150. (d) Critical value of the 
compression rate (a — a) /a as a function of k for A'' = 150. 



average body length L decreases with F monotonously. The snaking pattern appears for F > Fc, when the body 
length becomes sufficiently short. 

When an elastic rod is compressed, a buckling instability occurs if the compressive force F is larger than a critical 
value Fc- [l^jUll- In the classic Euler buckling of an elastic rod, a straight rod becomes unstable and a characteristic 
deformation of one- hump appears for F > Fc- The wavenumber of the deformation is q = n/L where L is the system 
size. The instability shown in Fig. 6(c) is interpreted as a kind of the buckling instability, however, a wriggle structure 
with many humps is characteristic in our system. 

We have numerically studied the wavenumber q of the snaking pattern in a larger system with N — 150. Figure 
7(a) shows a snapshot pattern for k = 47r/25 and F = 0.69. It is natural that the wriggle number of the pattern 
is larger in this larger system compared to the system of A'^ = 50 shown in Fig. 6(b). Figure 7(b) shows a similar 
snapshot pattern for k = 27r/25 and F — 0.69. The scale of the x and y coordinates is the same in Figs. 7(a) 
and (b). The wriggle number of the pattern is larger at k = 47r/25 than at fc = 27r/25. The average body size 
L = (a:jv/2 + a;Ar — a;i — /"^ is more contracted at fc = 47r/25 than at fc = 27r/25. We have studied the instability 

more in detail by changing the parameter fc. Figure 7(c) denotes the average wavenumber q of the snaking pattern 
as a function of fc for F — 0.4. Here, the average wavenumber q was numerically evaluated from the wavelength 
between the successive two peaks of the snaking patterns. The average wavenumber q is much larger than the value 
tt/Lq = 0.042 of the one- hump structure, where Lq = N /2 — 1. And, the average wavenumber q is definitely different 
from fc, but increases with the parameter fc. We have calculated the mean square {Sy^) and the average body length 
L = {xpf/2 + — xi — Xff/2+1) /"^ a-s a function of F at various fc's for N = 150 to find the critical points. The 
bifurcation occurs when the compression rate (a — a) /a is larger than a critical value. Here, a = L/{N/2 — 1) 
denotes the average interval between two neighboring elements in the a;-direction, and a = 1 is the natural interval. 
Figure 7(d) shows the numerically obtained compression rate at the bifurcation point as a function of fc for N — 150. 
The critical compression rate tends to increase with fc. We do not sufficiently understand why a snake pattern with 
q 2> tt/L appears in the ladder model. We neither understand the fc-dependence of q and the compression rate shown 
in Figs. 7(c) and (d). The analysis based on the theory of elasticity is left to future study. 



V. SUMMARY 



We have constructed several simple mechanical models which exhibit a uni-directional motion induced by internal 
forces and frictions. The model composed of two elements was first constructed for the bipedal motion. The energy is 
injected through the periodic change of the natural length of the spring, which plays a role of a kind of a muscle, and 
dissipated via the friction. This is a typical nonequilibrium system, although thermal fluctuations are not considered. 
The bipedal motion is caused by a kind of ratchet mechanism in which the reverse motion is suppressed by temporally 
enforced frictions. The sliding friction models might be more plausible, but the viscous friction models can be solved 
approximately. When the locomotive motion proceeds along an upward slope, the input energy is transformed to the 
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potential energy and a part of the input energy is dissipated by the friction. The efficiency at the maximum power is 
evaluated in these deterministic systems and it is compared to the theory of Van den Broeck. We have also proposed 
a modified model for the bipedal motion along a spatially periodic potential as a simple model of kinesin's motion. 
The effect of thermal noises was studied using the Langevin type equations. 

The bipedal models have been extended to a linear chain model and a ladder model. In the viscous friction model, 
we have found that Onsager's reciprocal relation is not always realized even if it is a linear system. In the ladder 
model, we have found a transition from the earthworm motion to the reptation motion when the forcing amplitude is 
increased. The transition occurs when the ladder is sufficiently compressed by the internal force and the friction. 

Our model systems arc very simple models to understand directional motions theoretically More realistic models 
including many degrees of freedom or three-dimensional deformations need to be developed to study specific systems 
such as molecular motors and real snakes. 
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